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C3 , ABSTRACT 



We present a detailed analysis of the X-ray emission of HS 1700+6416, a high redshift (z = 2.7348), luminous quasar, classified as 
a Narrow Absorption Line (NAL) quasar on the basis of its SDSS spectrum. The source has been observed 9 times by Chandra and 
once by XMM-Newton from 2000 to 2007. Long term variability is clearly detected, between the observations, in the 2-10 keV flux 
varying by a factor of three ( ~ 3 - 9 x 10~ 14 erg s _I cm -2 ) and in the amount of neutral absorption (N H < 10 22 cnT 2 in 2000 and 
2002 and N# = 4.4 + 1.2 x 10 22 cnT 2 in 2007). Most interestingly, one broad absorption feature is clearly detected at 10.3 + 0.7 
keV (rest frame) in the 2000 Chandra observation, while two similar features, at 8.9 ± 0.4 and at 12.5 ± 0.7 keV, are visible when 
the 8 contiguous Chandra observations of 2007 are stacked together. In the XMM-Newton observation of 2002, strongly affected 
by background flares, there is a hint for a similar feature at 8.0 ± 0.3 keV. We interpreted these features as absorption lines from a 
high velocity, highly ionized (i.e. Fe XXV, FeXXVI) outflowing gas. In this scenario, the outflow velocities inferred are in the range 
v= 0.12 - 0.59c. To reproduce the observed features, the gas must have high column density (N H > 3 x 10 23 cm -2 ), high ionization 
parameter (log£ > 3.3 erg cm s -1 ) and a large range of velocities (AV ~ 10 4 km s -1 ). This Absorption Line QSO is the fourth high-z 
quasar displaying X-ray signatures of variable, high velocity outflows, and among these, is the only one non-lensed. A rough estimate 
of the minimum kinetic energy carried by the wind of up to 18% L^/, based on a biconical geometry of the wind, implies that the 
amount of energy injected in the outflow environment is large enough to produce effective mechanical feedback. 
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r — , 1. Introduction als are observed in 10-15% of optically selected QSOs. Mini- 

00 ■ . . , , , , , BAL (500 <FWHM< 2000 km s 1 , Hall et al. 2002) and nar- 

^ . It is now widely accepted that most nearby bulge-dominated mw ab tion lines (NAL> FWHM< 500 km s 1 , Weymann et 
\Q : § a ™ m *i contain a ™ llc . Super Massive Black Holes aJ 19?9; G j gt aJ 2m) QSQ ^ observed in t0 30% 

1 A ' (S ^u ™LT T Tf 6nCe a tX f f tl0n Cen , m , aSS of optically selected QSOs. These absorption features are often 
in . of the SMBH and the properties (mass luminosity, ve ocity dis- j blue-shifted with respect to the source redshift, indi- 

^ i f r n slon) ° f * e h° st s P herold f ™ 6 "^ & Rlch St °f 199 lf eating outflowing velocities between 10 3 and few x 10 4 km s" 1 
£j , & Rees 1998, Magornan et al. 1998, Ferrarese & Merritt 2000), {Q J ^ & Bro f herton 2008) . 

. . ■ suggest that: 1) most bulge galaxies went through a phase of 

i> strong nuclear activity during their formation process and 2) the In the X-ray band, absorption due to ionized species such as 

' evolution of the galaxy and of the SMBH must be closely re- N VI- VII, O VII- VIII, Mg XI-XII, Al XII-XIII, Si XIII-XVI, as 

lated, implying some sort of interaction or feedback. Models well as L-shell transitions of Fe XVII-XXIV, is observed to be 

\ have shown that luminous AGN might sterilize their host galaxy, blueshifted by a few 10 2 - 10 3 km s _1 in ~ 50% of type 1 AGN 

by heating the interstellar matter through winds, shocks and ion- (the "warm absorber", e.g. Reynolds 1997; Piconcelli et al. 2005; 

izing radiation, strongly inhibiting star-formation in these galax- McKernan et al. 2007). In the last decade, thanks to the high 

ies and making their colors redder (Hopkins etal. 2008, Cattaneo collecting area of XMM-Newton, Chandra and Suzaku, blue- 

et al. 2009). Some observational evidence for such a scenario is shifted absorption lines due to highly ionized gas (i.e. Fe XXV, 

starting to be found in a handful of objects at high redshift (e.g. Fe XXVI) outflowing at high velocity (v „, = 1.5 - 6 x 10 4 km 

Cano-Diaz et al. 2012; Page et al. 2012). s~') have been observed in a number of local AGN (Markowitz 

One effective way of providing mechanical feedback can be et al. 2006; Braito et al. 2007; Cappi et al. 2009; Reeves et al. 

through massive, wide-angle AGN winds or outflows, that can 2009; Tombesi et al. 2010; Giustini et al. 2011). 

be observed as absorption features in optic al/UV/X-ray spec- , 

tra of AGN. In particular, UV broad absorption lines (BAL, , Flnall y' X " ra y BALs ' blueshifted up to v out ~ 2 x 10 5 km 

FWHM> 2000 km s 1 , Turnshek et al. 1980) of ionized met- s > have been observed in gravitationally lensed BAL QSOs at 

high redshifts (Chartas et al. 2002, but see also Hasinger et al. 

Send offprint requests to: Giorgio Lanzuisi, e-mail: 2002, Chartas et al. 2003, 2009, Saez & Chartas 201 1). These 

lanzuisi@iasfbo. inaf.it lines are associated with a zone of circum-nuclear gas, photo- 
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Table 1. Observation Log. (1) Telescope; (2) Detector; (3) Observation ID; (4) Observation date; (5) Exposure in ks (net exposure); 
(6) Net source counts in the 0.5-8 keV band. 



Instrument 


Detector 


ObsID 


Date 


Exposure (Net) 


Net Counts 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Chandra 


ACIS-I 


547 


2000-10-31 


50.1 (49.5) 


380 


XMM-Newton 


PN 


0107860301 


2002-05-31 


27.1 (10.8) 


312 


XMM-Newton 


MOS1 


0107860301 


2002-05-31 


27.1 (12.9) 


145 


XMM-Newton 


MOS2 


0107860301 


2002-05-31 


27.1 (13.2) 


155 


Chandra 


ACIS-I 


8032 


2007-11-12 


31.4 (31.0) 


205 


Chandra 


ACIS-I 


9757 


2007-11-13 


21.0 (20.7) 


124 


Chandra 


ACIS-I 


9756 


2007-11-14 


32.6 (32.3) 


221 


Chandra 


ACIS-I 


9758 


2007-11-16 


23.6 (23.4) 


194 


Chandra 


ACIS-I 


9759 


2007-11-17 


31.5 (31.2) 


201 


Chandra 


ACIS-I 


9760 


2007-11-19 


17.1 (16.9) 


153 


Chandra 


ACIS-I 


8033 


2007-11-20 


30.0 (29.7) 


204 


Chandra 


ACIS-I 


9767 


2007-11-21 


9.1 (9.0) 


60 
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Fig. 1. Chandra 0.5-8 keV, background subtracted, light-curve 
of HS 1700+6416 extracted from the 50 ks observation of 2000. 
The bin size is 3000 s. The dashed line shows the mean count 
rate value of 6.8 x 10~ 3 counts/s. 



ionized by the central X-ray source, with ionization parameter^ 
log£ ~ 3 - 5 erg cm s _1 and column density Nh~ 10 22 - 10 24 
cm -2 . These findings suggest the presence of massive clouds of 
highly ionized absorbing material, outflowing at nearly relativis- 
tic velocities from the nuclear regions, possibly connected with 
accretion disk winds and/or the base of a jet (e.g. Ghisellini et 
al. 2004). However, these features are difficult to observe, be- 
cause they are intrinsically variable, the gas is highly ionized 
(and thus only a few transitions are observable), and the spec- 
tral range where they occur lies at the high energy edge of the 
observing band available for present high throughput X-ray ob- 
servatories (Cappi et al. 2006). 

To understand the physical mechanism responsible for 
launching and accelerating AGN winds, it is crucial to build up 



1 £ = L/nR 2 , where L is the ionizig luminosity, n the particle density 
and R the distance between the ionizing source and the gas (Tarter et al. 
1969). 



a large sample of sources for which these X-ray features can be 
studied in detail, especially at high redshift, where the average 
SMBH accretion rate is higher (Barger et al. 2001). This will 
allow for the deeper understanding of the link between the ac- 
cretion and ejection processes, quantifying the kinetic energy 
injected in the environment, and thus the role of outflows in 
the formation of cosmic structures. In particular, good quality 
spectra of high-z QSO showing broad absorption features in the 
iron K band, are limited to three gravitationally lensed BAL 
and mini-BAL QSOs (APM 08279+5255, PG 1115+080 and 
H 1413+117, Chartas et al. 2002, 2003, 2007). These sources 
are the only few optical/UV BAL quasars to date that are bright 
enough in the X-ray band, mostly due to the strong lensing mag- 
nifications, to allow in-depth spectral analysis. 

Here we present a detailed analysis of the X-ray emission 
of the quasar HS 1700+6416, a high redshift (z=2.7348), high 
luminosity NAL QSO. HS 1700+6416 is one of the most lu- 
minous, non-lensed (Reimers et al. 1997), QSOs in the SDSS 
(Mi = -30.24). It is not detected in the NVSS radio survey, 
that has a 5 <x sensitivity limit of 2.5 mJy at 1.4 GHz. The 
source should be therefore considered extremely radio quiet 
(log(Fi. 4C //e /Fr < -5fl))- HS 1700+6416 is classified as a NAL 
QSO, showing narrow CII, CIV, Silll and Si IV absorption lines 
in the SDSS spectrum, blue shifted at mildly (~ 2.4 x 10 4 km 
s~') relativistic velocities (Misawa et al. 2007). The source has a 
very faint ROSAT counterpart (Reimers et al. 1995). 

The paper is organized as follows: Section 2 reports on the 
X-ray observations and data reduction; Section 3 presents the 
spectral analysis, while the main results are discussed in Section 
4. A standard A cold dark matter cosmology with Ho = 70 km 
s Mpc~', Qa = 0.73 and Q.m — 0.27 is assumed throughout 
the paper. Errors are given at 90% confidence level (AC = 2.706 
for one parameter of interest), unless otherwise specified. 

2. X-ray data analysis 

2. 1 . X-ray observations 

The quasar HS 1700+6416 lies in a particularly dense and inter- 
esting field in the sky. It is close to two clusters, Abell 2246, at 
z = 0.225 and 2.2' of angular separation, and V1701+6414, at 
z = 0.45 and 3.2' of angular separation (Vikhlinin et al. 1998), 



i z is the 1 .4 GHz radio flux density 
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and has a foreground over-density of galaxies and AGN, cen- 
tered at z = 2.300 + 0.015 (Digby-North et al 2010). For this rea- 
son, the field has been observed by Chandra and XMM-Newton 
several times (Tab. 1), with different scientific goals. In particu- 
lar, one 50 ks Chandra observation was performed in 2000, to 
study the cluster V1701+6414 (P.I.: L. Van Speybroeck); one 27 
ks XMM-Newton observation in 2002 was performed to study 
the same cluster (P.I.: F. Jansen), and a series of 8 Chandra ob- 
servations, for a total of ~ 200 ks, was taken during 9 days in 
2007, to study the rate of AGN activity among the galaxies in 
the over-density at z=2.300 (P.I.: K. Nandra). 

The 50 ks Chandra spectrum of HS 1700+6416 was pre- 
sented in Just et al. 2007, in a systematic study of the X-ray 
properties of the most luminous quasars in the SDSS. The au- 
thors applied the Cash statistic to the source unbinned data, and 
presented the spectra binned at a level of 10 counts per bin, 
for clearer presentation. In the spectrum of HS 1700+6416 they 
found two consecutive bins falling i 3<x below a simple power- 
law model, around 3 keV (11 keV rest frame). They claim the 
possibility of this being an absorption feature, and reported an 
improvement of the fit, with the addition of an absorption edge 
at 2.4 keV (9.0 keV in rest frame), significant at a confidence 
level > 99.6%. 

Misawa et al. 2008 presented an exploratory study of the X- 
ray properties of NAL-QSOs, analyzing the 50 ks Chandra spec- 
trum and the XMM-Newton spectrum. However they binned the 
Chandra spectrum to 20 counts per bin and the XMM-Newton pn 
and MOS to 150 and 45 counts per bin, respectively, at E> 2 keV. 
This led to the loss of any information about the presence of ab- 
sorption features in the spectra. Here we present the re-analysis 
of these two data sets, with additional 8 Chandra observations, 
with exposure times in the range 9-32 ks. 



2.2. X-ray data reduction 

The event files of the Chandra observations were retrieved from 
the Chandra X-ray Center (CXC) via the Web ChaSeR (Chandra 
Archive Search and Retrieval) interface. The data reduction was 
performed following the standard procedures outlined in the 
Science Analysis Threads for ACIS data at the CIAO (v 4.3) 
Web sit^B The 8 repeated observations of 2007 were meant to 
study the rate of AGN activity in the proto-cluster at z=2.3, and 
the quasar HS 1700+6416 was observed only as a by-product. 
As a result the source falls close to the ACIS-I ccd gap region, 
in all the 8 observations (with identical aim point). As a con- 
sequence, the effective total exposure at the source position is 
reduced by ~ 20% due to dithering. 

The XMM-Newton data were retrieved from the XMM- 
Newton archival Standard XMM-Newton S AS tasks epproc and 
emproc (SAS v. 7.1.00 were used to produce calibrated pn and 
MOS event files. The event files were processed using the latest 
calibration files and cleaned up from effects of hot pixels and 
cosmic rays contamination. X-ray events corresponding to pat- 
terns 0-12 (0-4) for the MOS(pn) cameras were selected. The 
XMM-Newton observation was almost completely affected by 
background flares (see the XMM-Newton User Handbook). We 
produced a global light-curve at energies >10 keV, where the 
contribution from the emission of real X-ray sources is min- 
imized. Different count rate cuts were performed, in order to 
find the best trade-off between increasing the number of source 




52000 52500 54415 54420 54425 
MJD MJD 




52000 52500 
MJD 



54415 54420 54425 
MJD 




52000 52500 
MJD 



54415 54420 54425 
MJD 



3 http://cxc.harvard.edu/ciao 

4 http://xmm.esac.esa.int 

5 http://xmm. esac. esa. int/sas/ 



Fig. 2. From top to bottom panel: Photon index T, column den- 
sity Neand 2-10 keV flux, in different observations. Time is in 
Modified Julian Date. The x axis scale of the right part of each 
panel is expanded to make visible the 8 contiguous observations 
of 2007. 



counts available and decreasing the signal to noise ratio. The fi- 
nal cut produced cleaned exposures of ~ 10 and ~ 13 ks for pn 
and MOS, respectively. 

We extracted the source counts from circular regions of ra- 
dius 4"for the Chandra pointings. We tried different extraction 
regions for the XMM-Newton EPIC cameras, in order to mini- 
mize the background contribution: 15", 25"and 40", i.e. 70, 80 
and 90% of the Encircled Energy Fraction (EEF) respectively. 
Given the relative flux level of source and background, we chose 
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the smallest extraction region (15"), that minimizes the contri- 
bution of the background, that could otherwise strongly affect 
the spectrum, especially above 2 keV. 

The background counts were extracted from annular regions 
close to the source, avoiding ccd gaps and other sources. The 
extraction areas were typically ~ 10 times the source extrac- 
tion region, in order to average out any position-dependent back- 
ground features. The obtained spectra have a number of source 
net counts, in the 0.5-8 keV band, of ~ 400 for the 50 ks 
Chandra observations, ~ 310(300) net counts for the pn (MOS1 
and MOS2 combined), and in the range 80-220 counts in each 
of the 8 repeated Chandra observations. 

Given the exposure times of the different observations, and 
the flux level of the source, the study of the short term variabil- 
ity is feasible only for the long Chandra exposure performed in 
2000. Fig. 1 shows the 0.5-8 keV background subtracted light- 
curve of HS 1700+6416, obtained from the 50 ks Chandra ex- 
posure, with a bin size of 3000 s, result in > 20 net source counts 
in each bin. When fitted with a constant, the resulting count rate 
is 6.8 x 10~ 3 counts s _1 , with^ 2 /ii.o./. = 21.4/16, and null hy- 
pothesis probability P(y 2 /v)= 0.16. Thus, the source is found 
consistent with being constant on time scales shorter than the 
duration of the observation (50 ks). 

3. Spectral analysis 

In the following, we will describe in detail the spectral analysis 
performed for the source: first, we determine the properties of 
the continuum, with a simple absorbed power law model; then 
we focus on the parameters of the strong absorption feature de- 
tected in at least 2 of the 4 available spectra, using a phenomeno- 
logical model (a Gaussian absorption line superimposed to the 
continuum); finally (Sec. 4) we investigate the possibility of the 
feature being produced by an absorption edge or by a highly ion- 
ized outflowing gas. 

Given the low number of source counts available in each ob- 
servation, and in order to investigate the presence and proper- 
ties of a possible absorption feature in the X-ray spectra of the 
NAL QSO HS 1700+6416, we performed all the spectral fits us- 
ing the Cash statistic (Cash 1979) implemented in Xspec vl2. 
This is a maximum likelihood function, that does not require 
counts binning, allowing the full exploitation of the Chandra 
and XMM-Newton cameras' energy resolution. The Cash statis- 
tic assumes that the error on the counts is purely Poissonian, and 
cannot therefore deal with background subtracted data. For this 
reason, a careful characterization of the background spectrum is 
needed, before performing a global fit to the source +background 
spectrum. 

The Chandra global background between 0.5 and 7 keV is 
complex, and has been reproduced with two power laws, mod- 
ified by photoelectric absorption, plus 3 narrow Gaussian lines 
to reproduce the features at 1.48, 1.74 and 2.16 keV (Al, Si and 
Au instrumental lines), plus a thermal (MEKAL) component, to 
account for the sum of the particle, cosmic and galactic com- 
ponents of the background (Markevitch et al. 2003, Fiore et al 
2012). The XMM-Newton EPIC pn (MOS) background between 
0.5 and 10 keV, has been reproduced with a similar model, but 
with 2 narrow Gaussian lines to account for the prominent fea- 
tures at 1.5 and 8 (1.8) keV (Nevalainen, Markevitch & Lumb 
2005). In each observation the background shape has been re- 
covered from background spectra extracted from a large (8 arc- 
min) region, excluding the astrophysical sources, in order to col- 
lect a large number of counts, thus allowing a detailed, global 
characterization of the background. Then, each model has been 



rescaled and fitted to the local background, leaving as free pa- 
rameters the overall normalization and the photoelectric absorp- 
tion, to account for position dependent variations. 

3. 1 . Continuum and Variability 

The source spectra have been fitted with a simple power-law, 
modified by intrinsic neutral absorption at the source redshift, 
plus galactic absorption (model wabs*zwabs*po in Xspec), fixed 
to the value measured by Kalberla et al. (2005) at the source 
coordinates0(jV H ,G fl /. = 2.3 x 10 20 ct«- 2 ). 

Fig. 2 shows the distribution of spectral parameters ifl NjJ^ 
and the 2-10 keV flux, as a function of time (MJD), for all the 
spectra available. The x axis scale of the right panel of each plot 
is expanded, with respect to the left panel, to make visible the re- 
sults obtained for each of the 8 contiguous observations in 2007. 

Long term 2-10 keV flux variability by a factor ~ 3 is clearly 
detected, from ~ 9 x 10~ 14 erg s _1 cirT 2 in 2002, to ~ 3.5 x 10~ 14 
erg s _1 crrT 2 in 2007. Also the amount of neutral absorption var- 
ied, from N H < 10 22 cnr 2 in 2000 and 2002, to 4 - 8 x 10 22 
cm -2 in 2007. The photon index is consistent with being con- 
stant, within the large error bars. 

The time span between the 8 observations performed in 2007 
is 9 days. In this period the spectral parameters of the source ap- 
pear to be constant, within the error-bars, allowing us to stack to- 
gether all the spectra in order to obtain an average, better quality 
spectrum. We summed all the source and background spectra, 
and produced weighted ARF and RMF with mathpha, addarf 
and addrmf tools, respectively, included in the HEASoft suiteQ 
The matrices are weighted on the basis of the exposure time, 
however we stress the fact that the 8 observations have exactly 
the same aim point, and thus the different ARFs and RMFs are 
almost equivalent. As a further check, we also stacked together 
the 8 event files, and extracted a single set of spectra and matri- 
ces, obtaining consistent results. 

HS 1700+6416 also has a rough mass estimate of Mbh=2.5x 
10 10 M (Shen et al. 2011). We stress that this extreme value, 
obtained from the C IV emission line, has large uncertainties, 
due to the possible presence of a non virialized (i.e. a disk 
wind) component in the CIV line (Shen et al 2008, 2011; Ho 
et al. 2012), and should be considered as an upper limit. We 
can roughly estimate the BH mass also from the 2-10 keV lu- 
minosity. The 2-10 keV intrinsic luminosity of the source goes 
from 7.27 x 10 45 erg s" 1 of 2002 to 3.85 x 10 45 erg s _1 of 2007. 
From the highest value, assuming that the BH is accreting at 
the Eddington limit, and using a bolometric correction factor of 
50 (Vasudevan & Fabian 2007), it is possible to derive a lower 
limit for the BH mass of 3.9 x 10 9 M that, together with the 
value reported in Shen et al. (2011) gives an indication of the 
possible BH mass range for the source. Considering the highest 
available mass estimate, we can derive a Schwarzschild radius 
(is = 2GMbh/c 2 ) of ~ 7.5 x 10 15 cm and a light crossing time 
of ~ 250 ks in the source frame, which correspond to a variabil- 
ity time scale of ~ 10 days in the observer frame, while for the 
lower value of the BH mass, the observed variability time scale 
is of 1.6 days. 



6 http://heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 

7 defined asFoc E~ l ', the power law spectral index. 

8 defined as M(E) = exp{-Ni,cr(E)), the equivalent hydrogen column 
density of the absorber, where <r(E) is the photo-electric cross-section. 

9 http://heasarc.nasa.gov/docs/software/lheasoft 
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Fig. 3. Left panel: Observed Chandra 50 ks spectrum and residuals of HS 1700+6416, modeled with a power law modified by 
Galactic absorption plus intrinsic neutral absorption. The counts are binned to a minimum significant detection of 3<x (and a max- 
imum of 15 counts per bin), only for plotting purposes. Right panel: Observed Chandra spectrum and residuals of the source 
obtained from the 8 observations of 2007. The same model and binning strategy of left panel are used. In red is shown the Chandra 
background. 



3.2. Absorption features 

Fig. 3, left panel, shows the spectrum of HS 1700+6416 obtained 
from the 50 ks Chandra observation. For plotting purposes only, 
the spectral counts are binned so as to have a minimum of signif- 
icant detection of 3<x, or 15 channel per bin otherwise. The fit to 
a simple absorbed power-law model shows significant residuals 
around ~ 2.7 keV (observed frame) suggesting the presence of a 
strong absorption feature. The C-stat value for this simple model 
is 137.3 for 170 d.o.f. When adding an absorption Gaussian line, 
with energy, line width and normalization as free parameters, 
the resulting C-stat is 117.0 for 167 d.o.f, i.e. a AC of 20.3 for 
3 additional parameters (we remember that the change in C-stat 
from one fit model to the next (AC) is distributed approximately 
as Ax 2 in this range of counts). 

Performing an F-test we find that the addition of an ab- 
sorption Gaussian line improved the fit at a confidence level 
> 99.999%. To further assess the statistical significance of this 
result, extensive (10000 trials) Monte-Carlo simulations have 
been carried out, producing simulated spectra with the FAKEIT 
routine in Xspec. The input model is an absorbed power-law 
with the best fit parameters obtained from the observed data. 
The simulated spectra are fitted at first with the simple power- 
law model, and then adding an absorption Gaussian line, with 
energy, width and intensity free to vary (range of Er, ne = 0.5 - 7 
keV). None of the simulated spectra shows a AC = 20 (the 
higher value being ~ 15) and therefore we estimate that the prob- 
ability to detect by chance the observed AC is < 1 x 10~ 5 (i.e. the 
feature is significant at confidence level > 99.999%, i.e. > 4<x), 
confirming the results obtained with the F-test. We stress the fact 
that, if the x 2 statistic and count binning (20 counts per bin) are 
applied, we obtain the following results: chi 2 /d.o.f '.=44.3/20 for 
the absorbed power-law model, and chi 2 /d.o.f .=17 .9/16, when 
an absorption Gaussian line is added. This turns into a slightly 
lower significance of 99.6%, in agreement with the value found 
in Just etal. 2007. 

As a further check we produced the ratio between the ob- 
served spectrum of HS 1700+6416, and the observed spectrum 



of a real featureless source, in order to exclude that some instru- 
mental effect could be responsible for the broad feature at 2-3 
keV. We selected the brightest source in the sample of BL Lacs 
observed in Donato et al. 2003, i.e PKS 2005-489, with a net 
count rate of 0.49 count/sec, because it has a featureless and flat 
spectrum and was observed with the same instruments (ACIS- 
I) and in the same period of the observation of HS 1700+6416. 
Fig. 5 shows the ratio between the two spectra. PKS 2005-489 
has a featureless, flat spectrum with F = 2.2, i.e. softer than 
the continuum of HS 1700+6416, as can be seen in Fig. 5 from 
the non flat ratio. More important, the feature between 2 and 3 
keV is clearly visible in the HS 1700+6416/PKS 2005-489 ratio, 
confirming that this is not an instrumental feature. 

The same analysis described above, was performed for the 
sum of the 8 Chandra observations of 2007. Fig. 3 right panel 
shows the stacked spectrum. The counts were binned as for the 
left panel, for plotting purposes. The fit to a simple power law 
plus neutral absorption model, in the 0.5-7 keV band, shows 
residuals around 2.4 keV, and a stronger feature at 3.3 keV 
(observed frame). The C-stat/d.o.f is 219.9/235 for the simple 
power-law model, while AC are 4.2 and 16.8 for the 2.4 and 3.3. 
keV features, respectively, for 3 additional parameters. This re- 
sults in a confidence level of 81.50% and 99.87%, respectively, 
when computed with F-test. 

We performed the same type of Monte-Carlo simulations 
adopted for the long Chandra exposure, finding that the prob- 
ability to detect by chance the observed AC is 1.03 x 10 1 and 
6 x 10~ 5 respectively (significance level 89.72% and 99.94%). 
Thus, the detection of the feature at higher energies is signifi- 
cant at > 3<r confidence level, while the presence of the feature 
at lower energies is only marginal (< 2a). 

Finally, Fig. 4 shows the EPIC pn and MOS spectra extracted 
from the XMM-Newton observation of 2002. The MOS1 and 
MOS2 spectra and responses were added with mathpha, addarf 
and addrmf tools. As already mentioned, the observation was al- 
most completely affected by background flares, and the cleaned 
exposure time is of ~ 10 ks. As a consequence, the number of 
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3.73 7.47 18.67 




Energy (keV) 

Fig. 4. Observed XMM-Newton pn (black) and MOS1+MOS2 
(red) spectra and residuals of the source, modeled with a power 
law modified by Galactic absorption plus intrinsic neutral ab- 
sorption. In green (blue) is shown the pn (MOS) background. 




- 0.5 1 2 5 

Energy (keV) 



Fig. 5. Spectrum of the ratio HS 1700+6416/ PKS 2005-489 for 
the 50ks Chandra observation (Cha2ooo in table 2). the spectra 
have been binned to ~ 200 eV bins. 



net source counts available in the 0.5-10 keV band is ~ 300 for 
each spectrum (pn and combined MOS). 

Despite the poor data quality, a hint of presence of an absorp- 
tion feature around ~ 2. 1 keV can be seen in the residuals of both 
pn and MOS cameras. In this case the C-stat/d.o.f is 493.2/529 
for the simple power-law model, while AC is 14.3 for 4 ad- 
ditional parameters (line energy, width and 2 normalizations). 
This results in a confidence level of 97.50% (F-test). Similar re- 
sults were obtained with Monte-Carlo simulations (98.15% con- 
fidence level). Thus the evidence of the absorption feature in the 
XMM-Newton spectra was marginal. 

4. Modeling of the absorption features 

We tried different models to reproduce the absorption features in 
the HS 1700+6416 spectra. The results are summarized in Tab. 



2 for the different models and data sets. Model 1 is the simple 
power-law, modified by intrinsic neutral absorption at the source 
redshift, plus galactic absorption, discussed above. 



4. 1 . Absorption Gaussian line 

Model 2 is the same model 1 plus one absorption Gaussian line 
(two for the Chandra 2007 spectrum). Even if this model rep- 
resent only a phenomenological model to assess the feature sig- 
nificance, the properties of the absorption lines, such as the line 
centroid energy, the width cr and the equivalent width (EW here- 
after), can be roughly recovered. 

The rest frame line centroid energy (Eu ne ) for the 50 ks 
Chandra spectrum is Eu„ e = 10.32^2 keV, the line width 
cr = 1.59+°-^ keV and the equivalent width EW = -0.83!^° 
keV (all quantities are rest frame). Fig. 6 left and right panels, 
show the 68/90/99% confidence contours of Eu„ e vs. cr and Eu ne 
vs. normalization, respectively, for the 50 ks Chandra spectrum 
(solid curves). Assuming that the absorption feature is due to 
either Fe XXV Ka or Fe XXVI Ka, the two strongest absorp- 
tion line for an highly ionized gas (Chartas et al. 2009, Saez et 
al. 2009), with rest frame energies of 6.70 keV and 6.97 keV 
respectively, we compute the minimum and maximum outflow- 
ing velocity, taking also into account the Icr error on E\- me . The 
resulting velocity is in the range v=0.30-0.46c. 

The observed width of the line in this spectrum is extremely 
large. This huge broadening can be due to either high turbulence 
velocity of the absorbing gas, or to the blending of a series of 
absorption lines, corresponding to different gas shells at differ- 
ent ionization states and/or moving at different velocities (or a 
combination of the three effects). We underline the fact that the 
energy of the Fe XXV Ka and Fe XXVI Ka lines being very 
close (A£'/, ne =0.27 keV), the simple blending of these two lines 
cannot be responsible for all of the broadening observed. 

The stacked Chandra spectrum of 2007 shows a weak fea- 
ture at ~ 2.4 keV and a more prominent one around ~ 3.3 keV. In 
the rest-frame, the line parameters are Eu ne \ = 8.86^43 keV and 
Eunei = 12.50+°;^ keV, oT < 0.95 keV and cr2 = 1.34+°;^ keV, 
EW1 = -0.14+°;* keV and £W2 = -0.5 lj° 2 f keV respectively. 
The dashed curves in fig. 6 represent the 68 and 90% confidence 
contours of the parameters of the two lines in the 2007 stacked 
Chandra spectrum. The resulting outflowing velocities are in the 
range vl= 0.20-0.3 lc and v2= 0.48-0. 59c, respectively. Finally, 
the line parameters for the XMM-Newton pn and MOS spec- 
tra fitted together are E Vme = 8.05+ ,- 3 , 3 , keV, cr < 0.42 keV and 
EW = -0.1 3^2 °4 keV. In this case the resulting outflowing ve- 
locity is v= 0.12-0. 22c. A hint of the presence of an emission 
line, close to the energy of the neutral Fe Ka line, is present 
both in the Chandra (2007) and in the XMM-Newton spectrum. 
However the feature is not statistically significant in either of 
them, and including such line in the spectral fit did not affect the 
properties of absorption feature. 

Fig. 6 clearly shows that the absorption feature is highly vari- 
able in energy, intensity and width. In particular, if we consider 
only the two features with the higher significance, i.e. the ones 
at 10.3 and 12.5 keV respectively, the probability that the en- 
ergy of the first absorption line is the same of the second is 
< 1 x 10~ 3 , i.e. the feature is variable in energy at a confidence 
level of > 99.9%. 
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Fig. 6. Left panel: 68,90,99% confidence contours of £/,•„<, vs. the line width from the APL+Gauss model fit of the Chandra 50 ks 
observation spectrum. The dashed (dotted) contours represent 68 and 90% confidence contours of the parameters of the two lines in 
the 2007 stacked Chandra (2002 XMM-Newtori) spectrum. Right panel: 68,90,99% confidence contours of En ne vs. normalization 
from the APL+Gauss model fit of the Chandra 50 ks observation spectrum. The dashed (dotted) contours represent 68 and 90% 
confidence contours of the parameters of the two lines in the 2007 stacked Chandra (2002 XMM-Newtori) spectrum. 



4.2. Absorption edge 

The third model we assumed to fit the spectral features of 
HS 1700+6416 was Model 1 modified by an absorption edge. 
The rest frame edge energy is E e d ge = 8.95^'^ keV and the 
maximum absorption depth is r = 1.85*2 g5 for the Chandra 
spectrum of 2000. For the stacked Chandra spectrum, the two 
edges have E edgel = 8.14^;| and E edge2 = U20$% keV, 
and rl = 0.60+°;^ and t2 = 0.82+°;^ respectively. Finally, 
the edge parameters for the XMM-Newton pn and MOS spectra 
are E edge = 7.63+°g keV an d t = 1.16+JJ ^. Given the differ- 
ent shape of the absorption edge, with respect to the Gaussian 
profile, the energies associated to the absorption edges are typ- 
ically lower than those found for a Gaussian line profile. The 
rest frame energies were always consistent with K shell ioniza- 
tion thresholds of mildly/highly ionized gas (Fe XVI-FeXXVI) 
with no blueshifted velocity (Hasinger et al. 2002), except for 
the highest energy edge of the stacked Chandra spectrum, at 
Eedgel ~ 11.2 ke V, that require an outflowing velocity of at least 
v/c=0.18 assuming the recombination energy of hydrogen-like 
iron, Fe XXVI of 9.28 keV. However, we don't expect to observe 
single edges in these kind of source, but rather complex absorp- 
tion structures, due to the chaotic nature of the gas, that should 
show a range of ionization parameters and velocities (Ebrero et 
al. 2011). 

4.3. Warm absorber 

As a further step we used the more physically motivated warm 
absorber model from XSTAR (Kallman & Bautista 2001), to 
modify the primary power law. The program computes the spec- 
trum transmitted through a spherical gas shell of constant den- 
sity, thickness AR, at a distance R >> AR from the ionizing 
source. In this analysis we used the program implementation that 
allows to include the XSTAR output directly within Xspec as a 
table/model. The input parameters of XSTAR are the gas tem- 
perature, density, luminosity and spectral shape of the ionizing 



continuum, turbulence velocity and element abundances, while 
the output Xspec model, superimposed to the absorbed power- 
law of Model 1, has 3 parameters that can be fitted to the spec- 
tral data: the ionization parameter the absorber column density 
Nh, and the absorber redshift (Model 4 in tab. 2). 

One important input parameter of the warm absorber model 
is the turbulence velocity of the gas Vturb> that controls the broad- 
ening of the absorption features. To reproduce the huge width of 
the absorption feature in the Chandra spectrum of 2000, we had 
to assume a large turbulence velocity, Vtwb = 3 x 10 4 km s _1 . 
This value is needed to mimic the observed spectral feature but, 
as said before, it can be interpreted as the result of different (con- 
tiguous?) gas shells moving with a large range of velocities, or, 
however, as the result of the probably complex dynamics of the 
wind. 

To reproduce the strong feature around the Fe line energies, 
without producing any detectable absorption at lower energies 
the absorber must have high column densities (N# > 5.5 x 10 23 
cm -2 ), and very high ionization parameter {log% - 3.85+i'H erg 
cm s _I ), such that almost all the Fe atoms are in the form of 
FeXXV-FeXXVI ions (Kallman et al. 2004), and most of lower- 
Z elements are completely ionized. In this model, the absorption 
feature is reproduced by a blend of the Fe XXV and FeXXVI 
Ka lines. Using the rest frame energies of these lines, i.e. 6.70 
and 6.97 keV respectively, the redshift of the absorber can be 
translated into a range of outflowing velocity comparable with 
that obtained from model 2 (v=0. 32-0. 45c). The improvement 
of the fit obtained with this model (AC ~ 19) is similar to that 
obtained with the Gaussian absorption line for the same number 
of free parameters. 

The same model was applied to the stacked Chandra spec- 
trum. The data require two highly ionized absorbers, with simi- 
lar column densities (~ 3 x 10 23 cm -2 ) and ionization parame- 
ters (Log(£) > 3.2 erg cm s _1 ) but outflowing at different veloc- 
ities, vl =0.19-0. 32c and v2=0. 47-0. 58c, to reproduce the fea- 
tures at 2.4 and 3.3 keV, respectively. In this case, the turbu- 
lence velocity required to reproduce the width of the lines is 
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Table 2. Results from spectral fits. Model APL is a power law with Galactic absorption and intrinsic absorp- 
tion (XSPEC model wabs*zwabs*pow); model APL+Gauss is an absorbed power law plus Gaussian absorption lines 
(XSPEC model wabs*zwabs(pow+zgauss)); APL+Edge is an absorbed power law plus absorption edges (XSPEC model 
wabs*zwabs*zedge*pow); model WARM(APL) is an absorbed power law modified by an intrinsic ionized absorbers (XSPEC 
model wabs*zwabs*warmabs*pow). All the energies are rest frame. Errors are given at 90% confidence level (AC = 2.706 for one 
parameter of interest, Avni 1976). 



Model 


Parameter 


Cha2ooo 


Cha2oo7 


XMM 


1 (APL) 


r 

N H (xl0 22 cm- 2 ) 
C-stat/d.o.f. 


1 98 +0 19 
< 1.20 
137.3/170 


9 ^q+0.22 

J -0.22 

4.43 + |' 2 ° 
270.1 /268 


2 06 +017 

^•^ -0.13 

< 0.64 
493.2/529 


2 (APL+Gauss) 


r 

N H (xl0 22 cm" 2 ) 


1.803H 
< 0.89 


2 - 20 -023 
2 95 + ™ 


2 04 +0 - 16 

-0 15 

< 0.41 




E tine (keV) 
o- (keV) 
EW (keV) 
C-stat/d.o.f. 


10 26 +0 72 
1 59 +077 

-0.83"S 
117.0/167 


8.86— 12.50— 
< 0.95 l-34 + « 9 

n , 4 +0.09 _Q 51 +0.23 
" -1 -0.12 "• J1 0.29 

247.0 262 


8 - 05 -a22 
< 0.42 

-0 13 +007 

478.9/525 


3 (APL+Edge) 


r 

N H (xl0 22 cm" 2 ) 


1 71+0.19 
i - ' 1 -0.19 

< 0.74 


1.98+0.25 

2 07 +i: ^ 


1 90+ 017 
< 0.35 




E edge (keV) 

T 

C-stat/d.o.f. 


8 95 +0 - 32 

i oc+0.56 

120.2/168 


a i A+o.n 11 9n+u.sv 

0.60+° 27 0.82*$ 
246.00 264 


7 63 +032 
-Q-2.2 

1 16 +m 

1 - lu -0.67 

482.8/528 


4 (WARM(APL)) 


r 

N H (xl0 22 cm" 2 ) 


1 89 +0 19 
< 1.03 


9 9 i +0.32 
< 3.63 


9 ac+0.31 

< 0.95 



z 1.50 ±0.15 1.87 ±0.14 0.94 ± 0.07 2.15 ±0.10 

N H (xl0 22 cm- 2 ) > 55.8 33.0^ 29.3^- 7 > 28.7 

log^Hergcms-') 3.85^ > 3.24 > 3.25 > 3.38 

C-stat/d.o.f. 118.3/167 247.6 262 482.7/527 



slightly lower, 5 x 10 3 and 10 x 10 3 km s -1 respectively. The 
goodness of the fit is again comparable to that obtained for the 
Gaussian absorption line model. When applied to the XMM 
spectrum, the model yields a column density Nh> 3 x 10 23 crrT 2 , 
Logf > 3.4 and v=0.12-0.20c, with an improvement of the fit of 
AC/Av = 10.5/3. 

4.4. Physical parameters of the gas 

Using the values of column density, ionization parameter and 
outflowing velocities derived from the XSTAR model for the 
different spectra, together with the estimated properties of the 
source, such as BH mass and X-ray (ionizing) and bolometric 
luminosity, we can give a rough estimate of the location, mass 
outflow rate and kinetic energy associated with the wind. We 
adopted the set of equations collected in Tombesi et al. (2012), 
and based on simple assumptions, like the definition of ioniza- 
tion parameter and escape velocity for the minimum and maxi- 
mum distance from the central source. 

The minimum radius, computed as the radius at which the 
outflowing velocity of the wind equals the escape velocity, varies 
in absolute values from the range 2.5 - 50 x 10 16 cm, assuming 
the very large (2.5 x 10 10 M ) BH mass estimate available from 
CIV, to 0.4-8.4x 10 16 cm, using the lower limit for the BH mass 
(3.9 x 10 9 M ), estimated assuming an Eddington rate accretion. 
In both cases (and independently of the BH mass) the very high 
outflowing velocities measured, imply an origin of the wind ex- 
tremely close to the BH, i.e. in the range 3-70 r s . The maximum 
radius can be inferred from the definition of ionization param- 



eter, as the maximum distance at which the ionizing luminos- 
ity can produce a ionization £ in a gas of column density N H : 
fmax = Lion/^Nn- The ionizing luminosity taken into account by 
XSTAR when computing the properties of the ionized gas is de- 
fined as the unabsorbed luminosity of the source in the range 
13.6 eV-13.6 keV. L,„„ is in the range 0.8 - 6.1 x 10 46 ergs s _ \ 
and the associated radius is comprised between 0.2 - 8.9 x 10 19 
cm, i.e. 2.6 x 10 2 - 1 x 10 4 r s for a BH mass of 2.5 x 10 10 M Q , 
and 1 .6 x 10 3 - 7.4 x 10 4 r s for a BH mass of 3.9 x 10 9 M . The 
values derived for the Chandra data of 2007 and the XMM data 
are upper limits, due to the lower limits in ionization parameter. 
The results obtained for r max are consistent with what found in 
Tombesi et al. (2012), for a sample of lower luminosity, local 
Seyferts, while the higher outflowing velocities observed in HS 
1700+6416, with respect to their results, translates into lower 
values of r min , i.e. the wind can originate closer to the BH. 

To compute the mass outflow rate we used the expression 
derived by Krongold et al. 2007 for a biconical wind, also used 
in Tombesi et al. 2012, M out = 0.%nm p NHV O utrf(5, <p). f{8, (/>) de- 
pends on the relative orientation of the disk, the wind and the 
line of sight, but for reasonable angles can be considered of the 
order of unity. The formula assumes a conical geometry with 
thickness constant and negligible with respect to the distance of 
the wind from the central source. The minimum mass outflow 
rates, associated with the minimum radius estimated above, are 
M out (r min ) = 4 - 6M Q /yr, and the associated kinetic energies 
are in the range E K {r mi ^) = 0.3 - 27 x 10 46 erg s" 1 . These val- 
ues translate into a minimum M ou t/M acc = 0.03 - 0.13 and a 
minimum E K /Lb i = 0.01 - 0.18, computed assuming for each 
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observation, the Lb i derived from the observed X-ray luminos- 
ity, an X-ray bolometric correction of 50, an efficiency 77 = 0.05, 
and constant density and velocity. Unfortunately, the maximum 
mass outflow rates remain unconstrained, due to the combina- 
tion of lower limits in Nh and However, such high minimum 
values are already very large, telling us that, if persistent and 
corresponding to the simple geometry assumed here, the mas- 
sive, fast outflowing wind observed in HS 1700+6416 carries 
a huge amount of energy, and is perfectly capable of provid- 
ing the mechanical power required to produce a significant feed- 
back in the environment, estimated to be of the order of ~ 5% 
of the bolometric luminosity (Di Matteo et al. 2005; Ostriker 
et al. 2010), or even a order of magnitude less in the case of 
multistage feedback (Hopkins & Elvis 2010). We also underline 
that the highest values of kinetic energy and mass outflow rate 
come from the feature with the most secure detection, i.e. from 
the Chandra spectrum of 2000, making our results more robust. 
Understanding the physics of these winds is necessary to obtain 
tighter constraints on their properties and on their impact on the 
surrounding environment. 



5. Conclusions 

We performed a uniform and comprehensive analysis of all the 
Chandra and XMM-Newton observations, taken between 2000 
and 2007, and determined the X-ray properties of the high red- 
shift NAL-QSO HS 1700+6416. The source was found to be 
constant in flux and spectral parameters on timescales of tens 
of ks, i.e. the duration of the longest observation available. 
Furthermore, the X-ray spectral parameters were constant during 
the 8 contiguous observations of 2007, performed in a time in- 
terval of 9 days. This allowed the stacking of the 8 spectra taken 
in 2007, resulting in a larger signal-to-noise spectrum for that 
period. The BH mass estimate available from the SDSS spec- 
trum is extremely high {M BH -2.5 x 10 10 M & ) and translates into 
a variability time scale of ~ 10 days in the observer frame. Using 
the bolometric luminosity and Eddington limit accretion, we de- 
rive a lower limit for the BH mass of 3.9 x 10 9 M e , implying a 
variability time scale of ~ 1.6 days. 

On longer timescales, the source is found to be variable on 
timescales of years, both in the observed 2-10 keV flux, that 
varies by a factor of 3 between 2002 and 2007, and in the amount 
of neutral absorption, that is negligible in 2000 and 2002, and 
becomes of the order of 5 x 10 22 ctrT 2 in 2007. 

The most remarkable feature however is the detection of a 
strong, broad and variable absorption feature, observed at en- 
ergies higher than that of the neutral Fe Ka emission line and 
up to ~ 12 keV. The strong, broad absorption feature is clearly 
detected, at high significance (> 4<x), in the spectrum extracted 
from the longest observation available, at an energy of ~ 10.2 
keV. In the stacked spectrum of 2007 a similar feature is detected 
(> 3<x) at higher energies (~ 12.5 keV). Furthermore, a hint of 
a feature at ~ 8 - 9 keV is observed, albeit at lower significance 
(~ 2cr), both in the stacked Chandra spectrum of 2007 and in the 
short (~ 10 ks) XMM-Newton spectrum of 2002. Finally, a clear 
variability is detected in the energy of the line, at least between 
the two strongest features observed at 10.2 and 12.5 keV (fig. 6). 

In addition to the phenomenological model (a power-law 
modified by a Gaussian absorption lines) used to determine the 
significance of the features, we also used a model including ab- 
sorption edges instead of Gaussian lines. The energies obtained 
are consistent with ionization thresholds of mildly/highly ion- 
ized gas (Fe XVI-FeXXVI) with no outflowing velocity, except 



for the highest energy edge of the stacked Chandra spectrum, 
for which a velocity v~ 0.2c is required. 

Finally, we adopted a more physical approach, modeling the 
absorption with the XSTAR photo-ionization code. The fit with 
XSTAR requires, for all the spectra, a dense (Nh> 3 x 10 23 
cm -2 ) highly ionized (log£ > 3.2 erg cm s _1 ) and highly tur- 
bulent (v filr i=0.01 -0.1c) absorber, such that a blend of Fe XXV 
and Fe XXVI Ka absorption lines reproduces the spectral fea- 
ture. The large turbulence velocity is required to reproduce the 
extreme width of the lines, however this broadening can be inter- 
preted as the result of the superimposition of different gas shells, 
with different velocities, rather than as an internal turbulence. 
The inferred outflowing velocities are in the range v=0. 12-0. 59c 

All these pieces of information make the NAL-QSO HS 
1700+6416 the best candidate to be the first non-lensed quasar 
showing X-ray BALs. The measured properties of this source in 
velocity range, line width and inferred physical parameters of 
the obscuring gas, are comparable only with those reported for 
the X-ray BALs APM 08279+5255 and PG 1 1 15+080 (Chartas 
et al. 2002, 2003). In particular APM 08279+5255 has BH mass 
(M BH ~ 2 x 10 10 M Q ) and Eddington ratio (L Bot /L Edd ~ 0.3) 
very similar to those of HS 1700+6416. Also similar is the large 
difference between the X-ray and UV absorption systems out- 
flowing velocities, (i.e. a ratio of the order of 10-15) possibly 
implying that UV and X-ray BAL are produced by different ab- 
sorbers at different distances from the central BH. 

With the parameters obtained from the fit comprising the ion- 
ized absorber, and assumptions on the geometry of the wind and 
the efficiency of the accretion, we were able to derive a rough es- 
timate of the location of the wind, that spans the range 3 — 10 4 r s , 
and of the minimum mass and kinetic energy content of the wind 
itself, that varies from few % up to ~ 13 and ~ 18% of the mass 
accretion rate and Bolometric luminosity respectively. Even if 
only approximated, these results, together with others, emerging 
from the analysis of local sources as well as from large blind 
search surveys, are starting to quantitatively probe that high ion- 
ized, massive, fast outflowing disk winds can provide the me- 
chanical energy needed to produce the feedback that might be 
responsible for the well known relation between galaxies prop- 
erties and BH masses. 
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